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ABSTRACT 

In this paper we present new empirical radio surface brightness-to-diameter (X — D) 
relations for supernova remnants (SNRs) in our Galaxy. We also present new theo- 
retical derivations of the T, — D relation based on equipartition or on constant ratio 
between cosmic rays and magnetic field energy. A new calibration sample of 60 Galac- 
tic SNRs with independently determined distances is created. Instead of (standard) 
vertical regression, used in previous papers, different fitting procedures are applied to 
the calibration sample in the log X — log D plane. Non-standard regressions are used to 
satisfy the requirement that values of parameters obtained from the fitting of X — D and 
D — X relations should be invariant within estimated uncertainties. We impose symme- 
try between X — D and D — X due to the existence of large scatter in both D and X. 
Using four fitting methods which treat X and D symmetrically, different X — D slopes (3 
are obtained for the calibration sample. Monte Carlo simulations verify that the slopes 
of the empirical X — D relation should be determined by using orthogonal regression, 
because of its good performance for data sets with severe scatter. The slope derived 
here (/3 = 4.8) is significantly steeper than those derived in previous studies. This new 
slope is closer to the updated theoretically predicted surface brightness-diameter slope 
in the radio range for the Sedov phase. We also analyze the empirical X — D relations 
for SNRs in the dense environment of molecular clouds and for SNRs evolving in lower- 
density interstellar medium. Applying the new empirical relation to estimate distances 
of Galactic SNRs results in a dramatically changed distance scale. 



Subject headings: ISM: supernova remnants — methods: analytical — methods: statis- 
tical — radio continuum: ISM 
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Introduction 



Supernova remnants (SNRs) are the main sources of kinetic energy and heat for the inter- 
stellar medium (ISM). They also contribute to the acceleration of cosmic rays. The total number 



1991 



of Galactic SNRs is p redicted to be between 1000 and 10,000 (jBerkhuiisenl 119841 ; iLi et al 
Tammann et al.lll994l ). Assuming a typical evolution timescale (~ 10 5 yr) of SNRs before they 



merge with the ISM an d an event rate of two supernovae (SNe) per century in the Milky Way 
(jDragicevich et al.lll999l ). ~ 2000 SNRs are expected in our Galaxy. However, just 274 SN Rs have 
been identified in our galaxy from their emission line spectra, radio and X-ray radiation ([Green 
20091 ). Such a great deficit points to the selection effects in the past surveys. All known SNRs are 
sources of radio-synchrotron emission. 

SNRs were often discovered in radio surveys of the Galactic plane. Because of the surface- 
brightness limit of previous surveys, more faint or confused SNRs await discovery. The dominant 
selection effects are those that are applicable at radio wav elengths. To put it simple, three selection 



effects apply to the identification of Galactic SNRs (e.g. lGreenlll99ll ): (i) difficulty in identifying 
faint remnants, (ii) difficulty in identifying small angular size remnants and (iii) the absence of 
uniform coverage of the sky. Additionally, Malmquist biaa^ acts severely in the Galactic samples 
making them incom plete. In the case of extragalactic SNRs, the selection effects are different (see 
Urosevic et al.ll2005l ). 



The important and very difficult task of determining distances to Galactic SNRs is often 
based on observations in radio domain. It is possible to determine direct distances from historical 
records of supernovae, proper motions and radial velocities, kinemat ic observatio ns, coincidences 
with H I, H II and molecular clouds, OB associations and pulsars (IGreenl Il984l ). Where direct 
distance determination is not possible, estimates can be made for shell SNRs using the radio 
surface brightness-diameter relationship (S — D). The mean surface brightness at a specific radio 
frequency, is a distance-independent parameter and is an important characteristic of an SNR 
(j Shklovskvl 1 1 960al ) . The radio surface brightness at frequency v to diameter relation for SNRs has 
been written in general form as: 

E„(Z>) = AD-?, (1) 

where A depends on the properties of the SN explosion and interstellar medium, such as the SN 
energy of explosion, the mass of the ejected matter, the density of the IS M, the magnetic field 
stren gth, etc., while (5 is thought to be independent of these properties (jArbutina h Urosevic 
20051 ). For empirical relations, parameters A and /3 are obtained by fitting the data for a sample 
of SNRs of known distances. 

According to the existing convention, SNRs are classified int o three basic types: shell rem - 
nants, plerionic or filled-center remnants, and composite remnants. Allakhverdivev et al. ( 1983al lbl. 



1986a 



b) showed that S — D relation is only applicable to shell- type SNRs. However it can also be 



1 The volume selection effect — brighter objects are favored in flux density limited surveys. 



-3 - 



used for composite remnants if we separate the flux from the surrounding shell and the flux from 
the central regions. 



The first theoretical E — D relation was derived by IShklovskvl (|l960al ) . H e was also the firs t 
to propose the use of this relation as a method for determining SNR distances (Shklovskv 1960b) 



Berezhko &; Volk 



An up dated theoretical derivation of this relation for s hell-like SNRs was done b y_ 
( 20041 ) . The first empiric a l relat i ons we r e deri ved by Poveda k, Woltier ( 1968 ). Milne! ( 197ol ) and 
later by IClark fc Caswelll (| 19761 ) . iMilnd ()1979l ) and many other authors. An important Galactic 
relation, frequently used in the last decade, was derived by C ase fc Bhattach a rya ( 1998, hereafter 
CB98 ). The updated Galactic E — D relation were derived bv lGuseinov et al.l (|2003l ) and lXu et al 



(|2005l ). More than five decades after the first published study of the relation, it continues to evolve 
both in theoretical and empirical aspects. A detai led review of Galactic and extragalactic empirical 
E — D relations was presented by lUrosevid (|2002l ). 



Although many authors have improved the theoretical and empirical E — D relations, th ere are 



still so me doubts when this relation is applied to determine distances to individual remnants. iGreen 
(j 19841 ) presented a critical analysis, noting that remnants with reliable distances have widespread 
intrinsic properties leading to significant dispersion in the T, — D diagram. One of the main problems 
is the high level of uncertainty for independently determined distances of SNRs, which are used 
for E — D calibration. Also, a variety of selection effec t s may be present in the data samples for 
both Galactic and extragalactic SNRs. lArbutina et al.l (120041 ) concluded that evolutionary paths 
may differ from remnant to remnant because of the potentially wide range of intrinsic properties 
of supernova (SN) explosions (and progen itor stars) and the interstellar medium (ISM) into which 
they expand. lArbutina Sz Urosevid (|2005l ) showed that different E — D relations can be constructed 
for two different samples, remnants in a dense environment of molecular clouds and SNRs evolving 
in ISM of lower density. 

A similar relation can be also used for composite SNRs, but only few of those remnants have 
known distances, and it is difficult to obtain a E — Z) relation for them. The E — D relation is 
often the only method for distance determination for shell SNRs. Despite inherent uncertainties 
and assumptions of this method, we suggest that the often-used E — D relation by CB98 should be 
updated because of the significant increase in the number of SNR calibrators. Additionally, different 
fitting procedures have to be used for establishing proper E — D calibration - we emphasize that 
in all previous papers the standard fitting procedure, based on vertical (parallel to y-axes) \ 2 
regression, was being used to derive empirical E — D relati ons for the Galactic SNR samples (e.g. 



Poveda fc Woltier! 1 19681 : IClark & Caswelll Il976l : lMilndll979l : CB98; iGuseinov et all 



2003 



Xu et al 



20051 ). In this paper, we select calibrators from Green's SNR catalog (IGreenl 120091 ) based on 
literature up to the end of 2008 H 



2 A Catalogue of Galactic Supernova Remnants (2009 March version), Astrophysics Group, Cavendish Laboratory, 
Cambridge, United Kingdom (also available on the World Wide Web at |http: / / www.mrao.cam.ac.uk / surveys /snrs / \ . 
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The £ — D relation has two particular uses: first, the derivation of diameters (and hence 
distances) and second, parametrization of the relationship between surface brightness and diameter, 
for comparison with models and theories. For the former, the measured value of £ is used to predict 
D, and for this case a least squares fit minimizing the deviations in log-D should be used. However, 
Urosevic et al.l (|2010l ) suggested the use of orthogonal regression fitti ng procedure for ob taining the 
empirical £ — D relation, instead of the classical vertical regression. Isobe et al. (jl99ol ) show that 
regression in astrophysical applications is more complex than most reali ze and give alterna t ives to 
ordina ry v ertical regression. In our paper, following the conclusions of llsobe et al .1 (|1990l ). iGreen 
(|2005l ) and lUrosevic et al.l (|2010l ). we test four different linear regression methods (including already 
mentioned orthogonal regression) that treat variables £ and D symmetrically. Symmetry between 
£ — D and D — £ is also required due to the existence of large scatter in both D and £ axes. 
Therefore, the values of parameters obtained by fitting £ — D and D — £ relations have to agree 
within estimated uncertainties. 

This paper also considers possible depende nce of radio luminosit y on linear diameter (L — D 
dependence). This criterion was established by lArbutina et al.l (120041 ). and is as follows: if L — D 
relation is obtained then £ — D relation follows and it may be used for the estimation of SNR 
distances. 

In Section 2, we present a new theoretical derivation of the £ — D relation based on equipar- 
tition between cosmic rays and magnetic field energ>H. The equations derived here describe three 
consecutive phases of the adiabatic (Sedov) phase of evolution for non-thermal radiation from SNRs. 
In Section 3, of the 274 SNRs in Green's present catalog (2009 March version), we construct a cal- 
ibration sample containing 60 Galactic shell remnants with distance estimates, and use it to derive 
£ — D relation. Additionally, Monte Carlo simulations are performed in Section 4 to estimate the 
influence of data scatter on the £ — D and L — D slopes. In Section 4, we also derive the empirical 
Galactic £ — D relation for the calibration sample containing 60 SNRs. In Section 5, we apply 
our Galactic £ — D relation to obtain the distances of two newly discovered, large and faint SNRs, 
G25. 1—2.3 and G178.2— 4.2. In Section 6, we present £ — D relations for a subsample of 28 Galactic 
SNRs evolving in dense environments and for a subsample only 5 Galactic SNRs evolving in low 
density. In Section 7, we briefly review selection effects which influence Galactic samples of SNRs 
severely. The conclusions of this paper are presented in the last section. 



2. A theoretical interpretation of the £ — D relation 

A slightly different theoretical interpretation of the T, — D relation is pre sented for non-thermal 
radiation from SNRs. This model is based on equipartition introduced by iRevnolds &: Chevalier 



3 By "equipartition" we do not assume that the energies are necessarily equal or nearly equal, but rather that the 
energy ratio is constant during evolution. 
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(|198ll ). The equipartition is a useful tool for estimating magnetic field strength and energy con- 
tained in the magnetic field and CR particles using only the radio synchrotron emission of a source. 
Details of equipartition and re vised equipartition calc ulatio ns for radio so u rces in general are 
available in iPacholczykl (|1970l ). iGovoni fc Ferettil (|2004l ) and iBeck &: Krausd (|2005l ). respectively. 
Arbutina et al.l (|2012r j introduced new modifications in equipartition calculation. They integrated 
over momentum to obtain energy densities of particles so there is no need to introduce a break in 
the energy spectrum and also take into account different ion species. 

The total energy density of cosmic rays (CRs) can be written as the sum of two components 
(j Arbutina et al.ll2012l ): energy of electrons (e e ) and ions (ei on ) 



ecR = e e + ei, 



kK c (m c 



2^2-7 



m p \ 2-7 
m. 



2m p c 2 E in j 
Ef n - + 2m e c 2 £ inj 



(7-l)/2 



r ( : V) r ( :i 2 1 ) \ - A (3- ) 2 , 

2V5F( 7 -l) ^ 1 



(2) 



where 7 is the energy spectral index (2 < 7 < 3), is injection energy of particles, m e and m p are 
electron and proton masses, V{ = rii/n are ion abundances, Ai and Z{ are mass and charge numbers 
of elements, and K e is the constant in the power-law energy distribution for electrons, k is a slowly 
varying function (if is not very high, see Fig. 2 in Arbutina et al. 2012), which incorporates 
the ratio between electron and ion energy (the latter of which should be dominant). Here, the term 
"total energy density" implies that we take into account all CR species (e.g. electrons, protons, 
a-particles, heavier ions) inside the SNR, which have been injected into the acceleration process. 
We neglect energy losses. 

We assume that an is otropic power- law distribution of ultra-relativistic electrons is created at 
the shock wave (jBelll Il978l ) 

N(E)dE = K c E~^dE. (3) 



For a given element of gas, K e evolves behind the shock wave as a result of conservatio n of energy. 



The fl ux density of synchrotron radiation of ultra-relativistic electrons, obtained from iPacholczvk 
(jl970l ) after substituting the emission coefficient e v with flux density S u , is 



S„ oc K c B l+a Vu 



W 



m 



2 Hz' 



(4) 



where B is the magnetic field strength, V is the volume, v is the frequency and a is the synchrotron 
spectral index defined as a = (7 - l)/2. Here, we use the flux density defined as 



(5) 
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where L v is the radio luminosity, / is the volume filling factor of radio emission, d is the distance 
and 9 = R/d is the angular radius. Then, we obtain the £ — D relation 



AD-? 



7T 



2 D 2 



oc S U D 



W 



m 2 Hz sr 



(6) 



where E„ represents surface brightness defined as = S u /Vt where f2 (in steradians) is the solid 
angle of the radio source. 



Reynolds fe Chevalierl ()198ll ) assumed redistribution of energy between cosmic rays, magnetic 



field and thermal gas behind shock wave (equipartition between ecR and e#) described with: 



ecR ~ = 7T— B 2 <*■ Po v s, 

ZflQ 



(7) 



where v s is the pre-shock velocity and po is the preshock ambient density. We introduce the param- 
eter s to allow for an ambient density profile of the form p oc R~ s , R oc t m an d 7 = 2a + 1, where 
R rep resents radius (D = 2R) and m is known as a deceleration para meter (Bandi era &; Petruk 
20041 ). This relation will also contain a break, similar to the relation of lDuric &: Seaquistl (|1986l ). 



If we assume low shock velocity (older SNRs) so that E[ n j 
7000 km/s, from equation ((2j), we have 



m P v s ^ 2m e c 2 i.e. v s <C 



K e oc p v s , 

and, after expressing B oc (/Oo^s) 1 ^ 2 from equation (|7|) and putting V oc D 3 , with 



(8) 



dt 



:0C 



mi 



m— 1 



oc mD~ 



we find the relation 



7+5 



S v OC K e 4 V" OC (po^s) 2 ^ 3 OC 



2(m-l) „ \ a + 3 



^+3 



0) 



Taking s = and m = 2/5 (|Sedovlll959l ) leads to 



3a + 7 

S„ oc L> — — , 



(10) 



and substituting an average S NR spectral index a = 0.5 gives j3 = 17/4 = 4.25 (see lRevnolds Sz Chevalier 



1981 



Berezhko fc VoUd 12004 ). 



Another possibility is E- m j 3> 2m c c 2 i.e. t> s 3> 7000 km/s, which is satisfied for younger 
remnants. In this case, following equation ([2]), constant K e in the power-law energy distribution 
for electrons has the form 
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K c OC po^ s ^inj « PO^, 



2(a+l) 



(11) 



and therefore the flux density is 



g + : J . 



3(a + l)(m-l) s(g + 3) 



+3 



(12) 



Taking s = and m = 2/5 yields 



9a + 7 

£„ oc D — ~, 



(13) 



and substituting a = 0.5 gives /3 = 23/4 = 5.75. 



Berezhko &; Volkl ([2004J) applied time-dependent nonlinear kinetic theory for cosmic ray (CR) 



acceleration in SNRs for studying the properties of the synchrotron emission. In particular, they 
applied detailed numerical calculations for deriving the surface brightness-diameter (Ti — D) relation 
for the range of the rel evant physical pa r amete rs, namely ambient density (njj) and supernova 
explosion energy (E sn ). iBerezhko fc Volkl (J2004J) derived (S — D) relations for different phases of 
SNR evolution. 

During the initial part of the free expansion phase, the expected dependence is 



(14) 



which explains well the numerical X„(D) behavior for small D. In this region Ti„(D) does not 
depend on the SN parameters E sn and M e j (mass of ejected matter), but does depend on the ISM 
density. 



During the later part of the free expansion phase the expected dependence is 



S, oc (£ sn /M cj ) 7/16 r4 1/16 ^ 5/16 



(15) 



In the subsequent Sedo\Q phase, we have 



E„ oc EU U D-^\ 



(16) 



independent of «h and M e j. 

In the late Sedov phase, dependence goes toward 



4 Also known as Sedov- Taylor phase 
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£„ oc E sn n^D- 2 . (17) 



We note that the phase described by equation (|13p . based on equipartition, represents an 
additional intermediate phase between late free expansion and early Sedov phase introduced by 



(jBerezhko fc Volkl 12004 )- Therefore, equations (fl~3j) . (fT6j) and (fTTl) respectively describes three 
consecutive phases of the adiabatic (Sedov) phase of evolution of SNRs. SNRs spend most of the 
time of their evolution in the Sedov phase, up to almost a million years in a hot ISM (McKee and 
Ostriker, 1977). Usually, pressure effects are responsible for terminating the Sedov phase earlier. 

Based on previous results and assuming that the majority of SNRs are in the Sedov phase, 
we expect to obtain £ — D slope j3 in the range of 2 to 5.75 depending on the evolutionary stage 
of the remnants in a sample of shell-type SNRs. Also, we should expect additional scatter (which 
does not affect the slope significantly) in the £ — D diagram due to different properties of the SN 
explosion and ISM, such as the energy of explosion, the mass of the ejected matter and the density 
of the ISM. 



3. Calibration sample and new £ — D relation 

To obtain a reliable empirical £ — D relation, we used a sample of shell-type SNRs with inde- 
pendently determined distances and well-determined angular diameters and flux densities. Various 
methods provide distance estimates to the remnants such as proper motions, shock and radial ve- 
locities, HI absorption and polarization, association or interaction with HI, HII and CO molecular 
clouds, X-ray observations, optical extinction and low frequency radio absorption. Of the 274 SNRs 
in Greens's present catalog (2009 March version), 84 have independently determined distances (59 
of them classified as shell type, 14 composite, 6 filled center and 5 unknown type). The flux den- 
sities given in the catalog and the calculated surface brightnesses in this paper are referenced to 1 
GHz. 

The surface brightnesses and angular diameters are taken from Green's catalog. This catalog 
was also our primary source of information about the distances of remnants. Many of the SNRs 
from this catalog have more than one distance available. Following CB98 approach, we have either 
chosen the most recent measurement or used an average of the available estimates (if the given 
distance range is narrow enough). We have also searched the literature for recent papers not 
included in Green's catalog, which provide accurate distances to Galactic shell remnants. Our 
Galactic sample of 60 shell remnants with direct distance estimates, which is used to derive new 
£ — D relation, is listed in Table 1. 

Part of our Galactic sample contains 33 SNRs, with updated parameters, which were also 
used as calibrators by CB98. The following 4 SNRs from the CB98 sample were omitted from our 
sample. According to Green's catalogue, the shell structure of G49.2-0.7 (W51) is questionable, 
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as well as flux density at 1 GHz, while G330. 0+15.0 (Lupus Loop) does not have enough accurate 
measurements of angular diameter and flux density. Up to now, only a lower limit for distance 
of rem nant G304. 6+0.1 (Kes 17) has been measured (9.7 kpc from HI absorption 



Foster et al 



(|2006l ) showed that OA 184, previously classified as SNR G166.2+2.5, is actually a Galactic HII 
region energized by 07.5V star BD+41°1144. 

Our Galactic sample also includes 5 composite SNRs because: (i) they have pure shell structure 
in radio regardless of centrally-brightened radio morphology; or (ii) it was possible to separate the 
shell flux density from the pulsar wind nebula (PWN) flux density. I Tarn et al.l (|2002l ) determined 
the integrated flux over various regions in the composite source Gil. 2-0. 3. From their measurements 
and integrated flux density measurements of Gil. 2-0. 3 at different frequencies, we estimate the SNR 
shell flux density at 1 GHz to be 20.5 Jy (subtracting PWN flux density which represents about 3 % 
of the flux density of the entire source). G93.3+6.9 (DA 530) is a high Galactic latitude SNR with 
a well-defined she ll-like radio morphology which has a centrally filled morphology only in X-rays 
(j Jiang et aJ 120071 ). Observed in radio, the composite SNR G189.1+3.0 (IC 443) consists of two 
connected, roughly spherical shells of radio synchrotron emiss ion, which are centered at different 
locations. Combining their measurements with existing data, ICastelletti et al.l (|2011al ) estimated 
that the flux density of the PWN represents only about 0.1% of total flux density and therefore we 
calculated the SNR shell flux density at 1 GHz to be 164.7 Jy. ICastelletti et al.l (|2011bl ) investigated 
in detail the radio emission belonging to SNR G338. 3-0.0 around the X-ray pulsar candidate in a 
search for traces of a PWN, by reprocessing data corresponding to observations acquired with the 
ATCA 0. No nebular radio emission has been found to correspond to the X-ray PWN, in spite 
of the good quality of their radio images down to low surface brightness limits, and therefore the 
measured flux density in radio corresponds only to the SNR shell flux densi t y. Ba sed on the radio 
images and the comparison with X-ray and IR observations, iGiacani et al.1 (120111 ) confirmed that 
there is no PWN within composite remnant G344. 7-0.1. Also, the same authors redetermined the 
distance of this SNR as (6.3 ± 0.1) kpc, on the basis of HI absorption and emission. 

The shell remnant Gl.9+0.3 is the youngest known Galactic SNR with known distance and 
the only Galactic SNR increasing in flux. Its estimated age is 156 ± 11 yr (jCarlton et al.ll201ll ) 
and distance about 8.5 kpc. Therefore, this SNR should not be included in our sample because it 
is probably still in the phase of free expansion. 

We mentioned in Section 1 that in the past a single linear regression method was used for the 
purpose of obtaining the empirical X — D relation: ordinary least-squares regression of the dependent 
variable Y against independent variable X, or OLS(Y|xfl In OLS(y|X), the regression line is 
defined to be that which minimizes the sum of the squares of the Y residuals. Some applications, 
however, require using alternatives to OLS(Y|X). The class of alternatives to OLS(y|X) used in 



5 The Australia Telescope Compact Array (ATCA), an array of six 22- m antennas used for radio astronomy 
6 Notation introduced by Isobe et al. (1990) 
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our paper was suggested by llsobe et al.l (|1990l ) for problems where the intrinsic scatter of data 
dominates any errors arising from the measurement process. This class of methods has also been 
proposed in order to avoid specifying "independent" and "dependent" variables. 

The most important purpose of the S — D relation is to estimate diamete r s (an d hence dis- 
tances) for Galactic SNRs based on their observed surface brightnesses. Green (j2005l ) also stated 
an important issue related to the £ — D fits, even in the case of disregarding problems with the 
selection effects. Namely, he pointed out that the measured value of £ is used to predict D, and 
that a least squares fit minimizing the deviations in log-D should be used to that effect. This, 
however, has not been done thus far, and fits minimizing the deviations in log £ have been used 
instead. He also mentioned that a fit to data that treats £ and D symmetrically is appropriate if 
a £ oc D n relation is used to describe the relationship between £ and D and, by extension, the 
radio-evolution of SNRs. 



Four methods that treat the variables symmetrically have been suggested bv llsobe et al.l (|1990l ). 
One is the line that bisects the OLS(Y|X) and the inverse OLS(X|Y) lines, called "OLS bisec- 
tor" or "double regression", which has been applied in characterization of the Tully-Fisher and 
Faber- Jackson relati ons to estimate galaxy distances ([Rubin et al.lll980l ; iLvnden-Bell et al.l Il988l : 
Pierce & Tullyill988l ). The second method is the geometric mean of OLS ( Y|X ) and the OLS(X\Y) 
slopes, proposed as the "impartial" regression by astronomer IStromberg ([l940h and used in cosmic 
distance scale applications. It was derived also by statisticians, independently, and called "reduced 
major-axis". The third regression is the line that minimizes the sum of the squares of the per- 
pendicular distances bet ween the data points and the line, often called "orthogonal regression" or 
"major-axis" regression. lUrosevic et al.l (|2010l ) applied orthogonal regression, for the first time, in 
obtaining empirical relation for SNRs in the starburst galaxy M82. The fourth method that treats 
the variables symmetr ically represents the arithmetic mean of the OLS(Y|X) and OLS(X|Y) slopes 
(lAaronson et al.lll986l ). It is easily recognized that these four techniques, though each is invariant 
to switching variables (referred to as the dependent and independent), lead to completely different 
regression lines, both mathematically and in real applications. 

The empirical £ — D relation depends greatly on the regression method adopted. The disper- 
sion of the six estimates is considerably larger than the variance of any one estimate (see Table 
2). For problems like these, Isobe et al. ( 199C ) suggest to calculate all sbi@ regressions and to be 
appropriately cautious regarding the confidence of the inferred conclusion. Accordingly, we give all 
mentioned regressions in Figure Q] (except OLS(X|Y), which is very close to orthogonal regression), 
applied to our calibration sample of 60 shell SNRs. 

Data fitting is performed analytically and by using bootstrap. The fit parameter values and 
their errors, presented in Table 2, were obtained after 10 6 bootstrap data re-samplings for each 



7 Actually, the sixth regression method, the arithmetic mean of the OLS(V|X) and OLS(X|F) slopes was not 
contained in their paper, but was added in proof. 
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fit. After applying bootstrap, we obtain the distribution of a selected value for a certain regres- 
sion procedure and then fit this dependence in OriginPro 8 with Gauss (normal) distribution, by 
using 4-parameter fitting without parameter fixing. A normal distribution is used here as a good 
approximation to describe the behavior of bootstrap re-samples. Thus, we get mean (expectation) 
and associated standard deviation for each parameter. 



4. Monte Carlo Simulations 
4.1. £ — D relation 

At first glance, inspection of Figure 1 and Table 2 leads to the conclusion that the resulting 
fit parameter values are significantly influenced by the type of fitting procedure. This is mainly 
due to the large dispersion in the sample of calibrators (low correlation coefficient^ r = —0.68). 
Much of this dispersion is due to properties of the SN explosion and ISM, which may substantially 
differ from one SNR to the other. It has generally been accepted that the density of the ISM is of 
significant importance in the e volution of SNRs; the other parameters are, in one way or another, 



connected to the ISM density (jArbutina fc Urosevid I2005I ) . SNRs of different types can be found 



along more or less parallel tracks in the £ — D plane. Also, the environment is probably quite 
inhomogeneous in the case of a single SNR, which would add more confusion in statistical studies. 
The errors of determined distances and, to a lesser extent, flux densities of the calibrator sample 
directly affect the scatter in the graph. 

Theoretical considerations assume that the T, — D relation corresponds to the evolutionary 
track of a typical SNR. Even for the Sedov phase, the fit should not necessarily be assumed to be 
linear (in log- log space) and Sedov sub-phases have (3 slopes from 2 to 5.75 as shown in Section 2. 
The previously described effects adds an extra scatter in our sample of calibrators. 

Due to the lack of information, it is not possible to separate all SNRs according to their intrinsic 
properties, which are connected with density of ISM in which they expand and their evolutionary 
stage. We therefore, expect to obtain, in some sense averaged, empirical £ — D relation that may 
still be used as distance estimator in cases when other existing methods are inapplicable. Errors of 
distances by this method will not be negligible but are still acceptable for the above purpose (the 
average fractional error of the CB98 relation defined as / = d ° h f~ dsd was about 40%). 

a obs 

We performed a set of Monte Carlo simulations to estimate the influence of scatter of the 
data on the T, — D slope. We generated random SNR populations (50 to 5000 SNRs) according to 
artificially chosen T, — D relations with /3 slopes 2, 3, 4 , 5 and 6, covering roughly the expected 
interval of slopes in Sedov phase. For each randomly selected value of D from interval 5 to 150 



Pearson product-moment correlation coefficient is calculated using the following equation: 
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pc (similar to that of real data) we calculate surface brightness £ using the initially proposed 
relation. Then we add random "polar" scatter in the created log-log data set. We add scatter by 
choosing random angle 4> an d random distance r from original point (x,y), after which we perform 
the following transformation (x, y) i — > (x + r cos cp,y + r sin (j)). For obtaining random angle (f> and 
distance r, respectively, we use uniform distribution and x 2 distribution with 10 degrees of freedom 
(usually labeled as xfo) implemented in numerical library for C programming language. 

We have chosen this distance distribution to be close to normal distribution, although additional 
tests showed that the type of distribution and number of degrees of freedom (for x 2 ) does not 
significantly affect our later conclusions. The only parameter which we enter is the distance that 
occurs most frequently in a generated distribution of distances (known as the mode in statistics) and 
therefore, it is a direct measure of dispersion from generated artificial samples (denoted by r mo( j e ). 
We have chosen such a mode value which ensures that the artificial distribution is similar to the 
real distributions for the 60 SNRs sample fitted with six different regression methods (r mo( je ~ 0.5). 

We give values of fitted slopes for randomly generated SNR populations in Figure i as a 
function of sample size and regression method (orthogonal, arithmetic mean, geometric mean, 
bisector, OLS(Y|X) and OLS(X|Y) regression). Each figure also contains information about the 
value of simulated slope f3 before performing random scatter on the log-log data set. 

After inspection of graphs on Figure [2] we come to the conclusion that orthogonal regression 
gives (3 slopes which are closest to the initial slopes of the artificially generated samples. Also, 
this type of regression shows significant stability as it converges to narrow value interval for larger 
sample sizes. 

This means that orthogonal regression is least sensitive to the scatter in the S — D plane. 
As mentioned at the beginning of this section, data dispersion is large in the calibration s e t tha t 



is used to construct our new empirical Galactic £ — D relation. Although llsobe et al.1 (jl990) 
concluded that the bisector performs significantly better than orthogonal regression, our class of 
problems (with large intrinsic scatter and steep slopes) gives more reliable solutions when solved 
with orthogonal regression. The existence of this significant scattering is obvious from the plots. 
It occurs as a result of coupling of several intrinsic SNR properties related to energy liberated by 
supernova explosions, density of surrounding media and evolutionary status of the remnants. The 
results of our Monte Carlo simulations (see Figure 2) show that the orthogonal regression give fitted 
slopes that are closest in value to the simulated slopes. Therefore, we strongly suggest orthogonal 
regression, instead of the previously used OLS(Y|X) and other regression types, for obtaining any 
kind of empirical £ — D relation with s lopes between 2 and 6. All previous Galactic empirical 



relations have slopes in this interval (see lUrosevid 120021 ) 



After applying non-weighted orthogonal regression on the sample containing 60 calibrators 
from Table 1, we obtain the relation 



'GNU Scientific Library (GSL), http://www.gnu.org/software/gsl/ 
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SigHz = 3.25™ • io- 14 ^- 4 -^ - 7 Wm- 2 Hz- 1 sr- 1 . (18) 

The new relation is significantly steeper than those obtained in earlier studies. The results of 
our Monte Carlo simulation (see Figure 2) also gave such a steep slope for the examined relation. 
Moreover, the orthogonal fit gives steeper slopes than any other regression, except OLS(X|Y). 
Shell SNRs with distances derived from our S — D relation are shown in Table 3. In total, 207 
remnants are shown. 

Next, we define fractional errors 



, (19) 
d\ 

in order to get an estimate of the accuracy of the £ — D relation for individual SNR distances 
and also as an indicator of the applicability of our relation for distance determination. Here d\ is 
the independently determined distance to an SNR and d-^ is the distance derived from the £ — D 
relation. Here, the average fractional error is / = 0.47 (comparable to that of CB98 / = 0.41, 
although we used a significantly larger sample). We also give average fractional errors for all six 
regressions in Table 2. 

The Monte-Carlo simulations presented here show that orthogonal regression is more stable 
for larger samples (see Figure 2), so our calibration data-set (60 SNRs) does not represent a proper 
sample for obtaining a high accuracy £ — D relation. A more stable behavior for larger samples is 
also the case for the other five regressions. 



4.2. L - D relation 

For a proper £ — D analysis, the L — D correlation should be checked. Using appropriate 
definitions of flux density and angular diameter, we have the following dependence £j, oc L U D~ 2 . 
However, the £ — D relation can be written as 



£ !/ = AD~ 2+d , (20) 

to allow for a possible dependence of luminosity on the linear diameter in the form L v = CD S , 
where C is constant. Radio luminosity as a function of the diameter for the calibrators in Table 1 
is given in Figure [3] and a very low correlation is evident. As an illustration, orthogonal regression, 
applied here, gives L — D slope 5 ~ —15.3. However, does not have a physical meaning because S 
should be around -3, according to equation (f20l) . Also, CB98 did not find any significant correlation 
between luminosity and linear diameter for shell type alone or for shell + composite-type distance 
calibrators. Previous studies of the £ — D relation indicate that 5 ~ 0, which leads to a radio 
luminosity that is independent of diameter and to the so-called trivial relation, S cx D~ 2 . If the 
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L — D correlation does not exist, the trivial E oc D 2 form should not be used ( Arbutina et al 



2004 ) 



Similar to the previous subsection, we also perform a set of Monte Carlo simulations to estimate 
the influence of the severe scatter of data on the L — D slope. We use a randomly generated E — D 
sample with added random polar scatter to construct L — D scatter by applying the luminosity- 
surface brightness relation L u = ir 2 D 2 T, u . 

Six regression methods were applied on randomly generated L — D samples. Figure H] shows 
fitted L — D slopes 5 as a function of sample size, with a scatter r mo d e = 0.5 and the simu- 
lated E — L> slope (5 = 5. At a first glance, Figure [4] leads to the conclusion that OLS(y|AT), 
bisector and geometric mean regression (average 5 values respectively —0.42,-1.41,-2.64) give 
much flatter slopes than arithmetic mean, orthogonal and OLS(X|Y) regressions (<5 respectively 
—8.84, —14.90, —17.25). We note that orthogonal regression, which is of significant importance for 
us, gives a very similar slope (—14.9) for randomly scattered data as that for real data (—15.3), 
which is very good agreement but with no physical importance. The orthogonal fit performs very 
unstably, especially for lower sample sizes (N < 500) and it is strongly influenced by data scatter, 
as are the other five fitting methods. Although OLS(Y|X), bisector and geometric mean regression 
seem more stable than other three methods, converting their slopes into angle interval leads to the 
opposite conclusion. 

The main reason for this instability can be inferred from brief inspection of Figure When 
we increase the scatter of the simulated data sets, the distribution of points tends to be almost 
vertical and even severe scatter can change the fitted slope from negative to positive. This effect 
produces a break in the dependence of the orthogonal fit slope from data dispersion (r mo d e ). This 
may be due to a lower span of luminosities for the L — D relation in comparison to the span of 
surface brightness for E — D relation (spans of diameters are the same for both relations). 

It can be inferred from further analysis that, taking into account all six fitting procedures, the 
geometric mean regression gives a slope that is close to that expected 5 = —f3 + 2 (for simulated 
T, — D slope (3 = 6,5,4,3). Figure [6] also shows that this type of regression can be a valid tool 
for obtaining L — D relation for a limited range of data dispersion. Thus, we should be careful if 
geometric mean regression is applied on L — D analysis since it is not insensitive to severe scatter 
like orthogonal regression when applied on E — D relation. 

Therefore, we suggest to obtain the L — D relation directly by using some type of regression 
only if the data set of calibrators is not subject to severe scatter. In other cases, a more accurate 
approach would first require the derivation of E — D relation from which then, by usi ng the equation 



201) . we obtain the L — D relation. Our conclusion is also supported by the fact that lArbutina et al 



(|2004l ) did not find any significant correlation when analyzing L — D dependence, except for the 



M82 starburst galaxy (with a fit quality of 64 per cent) for which a good E — D relation does exist. 
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5. Distances of two newly found Galactic SNRs: G25.1-2.3 and G178. 2-4.2 

The Sino-German A6 cm Galactic plane survey is a sensitive survey with the potential to 
detect new low-surface-brightness SNRs. This survey searches for new shell-like objects in the A6 
cm survey maps and studies their radio emission, polarization and spectra using the A6 cm maps 
together with the All cm and A21 cm Effelsberg observations. 



Gao et al.l ([201 ll ) have discovered two new, large, faint SNRs, G25. 1-2.3 and G178.2-4.2, 
both of which show shell structure. G25. 1—2.3 is revealed by its strong southern shell, which has 
a size of 80' x 30'. It has a non-thermal radio spectrum with a spectral index of a = —0.49 ± 0.13. 
G178.2— 4.2 has a size of 72' x 62' with strongly polarized emission being detected along its northern 
shell. The spectrum of G178.2— 4.2 is also non-thermal, with an integrated spectral index of a = 
—0.48 ± 0.13. Its low surface brightness makes G178.2— 4.2 the second faintest known Galactic 
SNR. This demonstrates that more large and faint SNRs exist, but are very difficult to detect. 



Using the total intensity and the radio spectral index, iGao et al.l ([201 ll ) calculated the surface 
brightness of both objects at 1 GHz, obtaining Eighz = 5.0 • 10 -22 Wm _2 Hz _1 sr _1 for the southern 
shell of G25.1-2.3 and Si G h z = 7.2 • 10" 23 Wm^Hz^sr- 1 for G178.2-4.2. 

The S — D relation of CB98 was used to estimate the distances of the two newly found SNRs. 
For G25. 1—2.3, they found that the diameter for the shell is 72 pc and the distance is 3.1 kpc, 
consistent with those derived from the HI data. For G178.2— 4.2, they found that the diameter 
is 197 pc and its distance is 9.4 kpc. This places the object far outside the Galaxy, which seems 
hardly possible. 

This is the appropriate place to mention specific behavior of the £ — D relation, depending 
on the brightness of the SNR for which the distance is calculated. The slope of the £ — D relation 
directly affects derived diameters and hence distances. A steeper slope will give larger diameters 
(and hence distances) for SNRs with higher surface brightnesses, while it will lead to smaller 
diameters (distances) for SNRs with lower surface brightnesses. 

We applied our new Galactic £ — D relation to obtain the distances of these two objects. For 
G25. 1—2.3, we obtained that the diameter for the shell is 42 pc and hence the distance is 1.8 kpc, 
which is still not far from that derived from the HI data. For G178.2— 4.2, we found that the 
corresponding diameter is only 63 pc, so that its distance is 3.0 kpc, which places the object inside 
the Milky Way. Note, however, that the uncertainties of these estimates could be as large as 50%. 



6. Dependence of £ — D relation on the density of the interstellar medium 



We focus on a mo re direct connection between the ISM density and the £ — D relation. 



Duric h Seaquistl (j 19861 ) assumed that the dependence of surface brightness on the density of the 



ISM has the form £ oc pi oc where po an d bh are the average ambient density and hydrogen 
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number density, respectively, and rj is a constant (see also iBerezhko Volkl 120041 ) . This means, 
the larger the surrounding ISM density, the greater the synchrotron emission from the SNR. 

This suggests that, on average, SNRs in dense environments would tend to have higher surface 
brightness in comparison to those evolving in lower-density medium. However, since X = AD~ 17 ^ 
relation depends primarily on the initial explosion energy (equation 1 16f) . variation in SN energy 
rather than density could produce scattering in some parts of the X — D diagram. 

Although it is difficult to decide which SNRs should be included or excluded from analysis, we 
tried to extract subsamples of Galactic SNRs (using Table 1) based on the density of the medium 
in which they evolve and then to obtain particular X — D relation. 



Mathweson et al.l (|1983l ) classified SNRs into four categories based on their optical characteris- 
tics: Balmer-dominated, oxygen-rich, plerion/composite and evolved SNRs. Balmer-dominated line 
emission is produced when the expanding SNR encounters low-density (~ 0.1 — 1cm -3 ), partially 



:< 10 



15 



cm lomza- 



neutral ISM. The hydrogen overrun by the shock is collisionally excited in a thin 
tion zone, producing optical spectra dominated by the Balmer lines of hydrogen (jGhavamian et al 



20071 ) . Oxygen-rich SNRs were mainly identified based on their optical oxygen line emission prop- 
erties and they primarily occur in HII and molecular cloud regions, i.e. in higher density regions. 
Young oxygen-rich SNRs are likely to interact with especially complex circum-stellar medium, 
rather than ISM. Balmer-dominated SNRs are connected to Type la SNe - deflagration of C/O 
white dwarf, while oxygen-rich SNRs originate in the Type lb events - explosions of a massive O 
or a Wolf-Rayet (W-R) star. 

We created a subsample of 28 Ga lactic SNRs with known distances, using online database as 



part of pape n I iGilles &: Samarl (12012 ). interacting with molecular clouds, i.e. evolving in dense 



ISM. The properties of these SNRs are given in Table 4. After applying orthogonal regression on 
the subsample from Table 4, the X — D relation obtained is 



X 



1GHz 



3.89™ 1 • W -^ D -3.9±0A Wm -2 Hz 



(21) 



which is slightly steeper than relation obtained by lArbutina et all (120041 ) /3 = 3.5 ± 0.5, whose 
sample contained only 14 Galactic SNRs associated with molecular clouds. The average fractional 
error is / = 0.35. Slope j3 = 3.9 for Galac tic SNRs in dense IS M is in very good agreement with 
slope obtained for the M82 data sample by lllrosevic et all (120101 ). w ho also used orthogonal fi tting. 
Our obtained slope approximately coincides with the conclusion of IBerezhko &: Volkl (|2004l ) that 
SNRs in dense ISM should populate the part of the line X = AD~ 17 / A . 

The sample of Galactic SNRs evolving in low-density ISM is significantly smaller. Up to now, 
only 5 Galactic SNRs showin g Balmer-dominated characteristics have been reported (see review 
papers IVinkll2012l : iHengJ 12010 ) and are therefore thought to expand in low-density medium. The 



10 This paper builds on the List of Galactic SNRs Interacting with Molecular Clouds maintained by Bing Jiang 
"http : //astronomy . nju . edu . cn/~ygchen/ others/bj iang/ inter SNR6 . htm" 
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Balmer-dominated sub-sample contains: Kepler (G4.5+6.8), Cygnus Loop (G74. 0—8.5), Tycho 
(G120.1+1.4), RCW 86 (G315.4-2.3) and SN 1006 (G327.6+14.6), whose basic properties can be 
found in Table 1. The £ — D relation obtained is 

SigHz = 1.89^29 • 10- 16 D- 3 - 5±0 - 5 Wm-^sr" 1 . (22) 

with an average fractional error / = 0.18. Figure \7\ shows the relation for Balmer-dominated and 
oxygen-rich Galactic SNRs for comparison purposes. The two relations have similar slopes, but the 
former is below the latter, as we would expect for SNRs in low-density environments. To further 
investigate the influence of the ISM density on the £ — D relation, we should use more complete 
samples, especially sample of Balmer-dominated remnants. 



Selection effects 



Identification of Galactic SNRs is always accompanied by selection effects ar ising from di f ficulty 



in identifying (i) fa i nt SN Rs and (ii) small angular size SNRs (e.g. CB98, iGreen 



Urosevic et al 



1991 



2005 



2005 



201CJ). A high enough surface brightness for SNRs is required in order to 



distinguish them from the background Galactic emission. Any Galactic radio survey is severely 
biased by observational constraints. The surface brightness limit seriously affects the completeness 
of all catalogues of Galactic SNRs. The derived flux densities of many SNRs are also poorly 
determined. Small angular size SNRs are also likely to be missing from current catalogues and 
calibration samples. Their small size may cause wrong classification of SNR type and also lead to 
incorrect calculation of their surface brightness. 

Additionally, Malmquist bias severely acts in the Galactic samples making them incomplete. 
This is a type of volume selection effect which naturally favors bright objects in any flux-limited 
survey because they are sampled from a larger spatial volume. The result is a bias against low 
surface-brightness remnants such as highly evolved old SNRs. Only the extragalactic samples are 
not influenced by the Malmquist bias since all SNRs in the sample are at the essentially same 
distance. 



Urosevic et al.1 (|2010i ) performed a set of Monte Carlo simulations to estimate the influence 
of the survey sensitivity selection effect on the S — D slope for the M82 sample (31 SNRs). An 
appropriate sensitivity cutoff is applied to the simulated data points, taking into account only 
points above the given sensitivity line. They concluded that the sensitivity selection effect does not 
have a major impact on the S — D slope. 

Technological advances in radio telescopes and X-ray instruments will greatly increase the 
number of known supernova remnants (SNRs), lead to a better determination of their properties 
and, thus, reduce the influence of the mentioned selection effects. 
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8. Summary and conclusion 

This paper presents a re-analysis of the theoretical and especially empirical Galactic £ — D 
relation and also the dependence of this relation on the density of the ISM. Motivation was found in 
the observed property that the empirical X — D relation strongly depends on the used regression due 
to severe scatter in data samples. In contrast to the standard least-squares (vertical) regression, 
we examine the behavior of six different types of regression. We put emphasis on those which 
satisfy the requirement that the values of parameters obtained from the fitting of £ — D and D — £ 
relations should be invariant within estimated uncertainties i.e. treat £ and D symmetrically 
(namely: orthogonal regression, bisector, arithmetic and geometric mean regressions). 

The catalog of known Galactic SNRs has grown in size significantly since the work of CB98, 
from 215 to 274 SNRs. The number of SNRs with known distances has also increased. We included 
the latest distance updates, if available, to derive a new Galactic £ — D relation, using a sample 
of 60 shell SNRs. We concluded from our tests that the size of the sample significantly influences 
the stability of any type of regression, and that more observations are necessary to increase the 
statistical significance of the sample. 

We present an additional modification to the theoretical £ — D relation for SNRs in the 
adiabatic expansion (Sedov) phase. This modification is based on the equipartition introduced by 



Reynolds &: Chevalieii (|198ll ). Our modification induces a new sub-phase lying between late free 



expansion and early Sedov phase having slope = 5.75. Also, ou r theoretical slope (0 = 4 .25) for 



relatively older SNRs in Sedov phase is in agreement with that of lBerezhko &: Volkl (120041 ) 



We have performed an extensive series of Monte Carlo simulations to evaluate numerically how 
well each of the six regression methods approximate the artificially proposed £ — D dependence 
with added scatter. The orthogonal regression is the most accurate slope predictor in data sets with 
severe scatter as it gives the slope closest to the original value. For other regression procedures, the 
fitted slopes of the £ — D relation are seriously affected by data scatter. Standard (vertical) fitting 
or OLS(y|X), which was used in previous papers considering empirical £ — D relation, leads to a 
significant change in slope when applied on highly dispersed data sets. When the scatter increases, 
the standard fitting leads to flatter slopes close to trivial (0 « 2). Also, similar conclusions apply 
to OLS(A|y), bisector, arithmetic and geometric mean regression. A possible reason may be the 
inherent property of £ — D data sets which are defined by a narrow range of diameters (about 1 
order of magnitude) in comparison to the wide range of surface brightness (approximately 5 orders 
of magnitude). We recommend using orthogonal regression for obtaining any type of empirical 
£ — D relation. 

Our obtained slope (0 = 4.8) is significantly steeper than previous values, it is far from the 
trivial one {0 ~ 2.0) and agrees with theoretical predictions for the Sedov phase of SNR evolution. 
Since our calibration sample contains 60 SNRs, which probably have different explosion energies, 
evolve in different ambient density and may be in different phases of SNR evolution, the relation 
£ oc D~ 4 8 could represent an averaged evolutionary track for Galactic SNRs. As such, it could be 
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used potentially for estimating the distances of SNRs in our Galaxy. A relatively large data scatter 
can be partly explained by the above mentioned influences. 

Additional analysis is carried to examine possible dependence of radio luminosity on linear 
diameter (L — D relation) for SNRs. Also, this analysis may provide an answer as to whether 
determination of SNR distances on the basis of a X — D relation is possible. We did not find any 
significant correlation for our sample consisting of 60 SNRs, while orthogonal regression gives a 
very steep slope (/? ~ —15) which has no physical meaning. Monte Carlo simulations have revealed 
that the L — D relation is very sensitive to severe scatter in the data, which is surely present in 
our sample. Therefore, it is possible to obtain L — D relation only if the set of calibrators is not 
subject to severe scatter. Otherwise, we obtain very steep slopes without physical meaning. We 
recommend an indirect approach to obtaining an L — D relation, which requires obtaining the X — D 
relation first (by orthogonal regression). 

We made an attempt to find a homogeneous subsamples of Galactic SNRs in low-density and 
dense environments and to obtain the X — D relations for two particular classes of SNRs. Applying 
orthogonal regression to the sample (which contains 28 SNRs) made of SNRs which evolve in dense 
environment of molecular clouds (including oxy gen-rich SNRs), l eads to the slope /3 = 3.9. Our 



sample contains twice as many SNRs as that of lArbutina et al.l (120041 ) and the obtained slope is 



sli ghtly steeper. Our slo pe for SNRs in dense ISM is in good agreement with the result obtained 



bvlUrosevic et al.l (|2010l ) for the sample of 31 SNRs in the M82 galaxy (/3 = 3.9). lUrosevic et al 
proposed their relation to be used for estimating distances to SNRs that evolve in a denser 



interstellar environment, with number density up to 1000 particles per cm 3 . This agreement may 
be very useful because this is the best sample for the X — D analysis, consisting of a relatively high 
number of very small and very bright SNRs from starburst galaxy M82. 

In contrast to the sample evolving in dense environment, we still have a very small sample 
of SNRs in low-density medium (Balmer-dominated), consisting of 5 SNRs only. Obviously, this 
sample is too small for any firm conclusions to be made. The X — D slope (3 = 3.5 obtained for this 
sample is close to that for SNRs in a dense environment and these two classes of SNRs approximately 
lie on two parallel tracks or domains in the log-log plane, one above another, as expected. However, 
the small number of objects in the Balmer-dominated sample strongly constrain the reliability of 
this relation, therefore it should be used with caution. More optical and X-ray observations are 
needed for discovering new Balmer-dominated (la) and oxygen-rich (lb) SNRs. 

We used our new empirical relation to estimate distances to 147 shell-like remnants with 
unknown distances and obtained a drastically changed the distance scale for Galactic SNRs. Though 
we are aware of theoretical and statistical flaws in the X — D relation, we think that in cases where 
direct distance estimates are unavailable, the T, — D relation should remain an important tool 
for distance determination. The obtained X — D relations should be used with caution, because 
uncertainties of distance estimates could be as large as about 50%. Although we have increased 
the number of calibrators significantly, more observations are needed to get a better understanding 
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of the radio evolution of SNRs. 
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Fig. 1. — The surface brightness vs. diameter £ — D relation at 1 GHz for shell SNRs obtained 
by using the distance calibrators in Table 1. The different methods for minimizing the distance 
of the data from a fitted line are presented. Two solid lines represent OLS(Y|X) (thin line) and 
orthogonal regression (thick line). Dash, dot and dash-dot lines represent arithmetic, geometric 
mean of the OLS(y|X) and OLS(X|Y) slopes and OLS bisector respectively. OLS(X|Y) line, with 
slope very similar to orthogonal regression, is omitted to avoid complicated graph. G1.9 + 0-3 , 
the youngest Galactic SNR which shows the flux density increasing with time (IGreen et al.ll2008l ). 
is also shown (triangle) but not included in th e calibration sample be cause it is still in the early 
(raising) free expansion phase of evolution (see iBerezhko &: Voildl2004l ). 
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Fig. 2. — The fitted S — D slopes vs. sample sizes for six regression types applied on randomly 
generated SNR samples containing artificial scatter and simulated slopes /3=2, 3, 4, 5, 6. This 
analysis was done to determine which type of regression is the least sensitive to scattering. 



1E14 I — — — — <-•■ 

1 10 100 1000 

D[pc] 



Fig. 3. — The radio luminosity (at 1 GHz) vs. the diameter for the SNRs ( £ — D calibrators) given 
in Table 1, to allow a possible dependence of the luminosity on linear diameter in form L v = CD 5 . 
Weak correlation is evident, primarily due to severe data scatter. 
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Fig. 4. — The fitted slopes of the power-law fit L u = CD vs. sample sizes for six regression types 
applied on randomly generated samples (r mo d e = 0.5). It is evident that the orthogonal fit becomes 
very unstable for L — D relation and strongly influenced by data scatter, especially for a small 
number of points, as is the case with our sample. 



-29- 



-14 



-15 



-16 



-17 



60 

o -18 



-19 



-20 



-21 



-22 



mode 



• 



• •• 



_l I I I I I I I I I I I I L 

-2 -1 1 2 3 4 5 



logD 



-14- 



-15 



-16 



-17- 



g> -18- 



-19- 



-20 



-21 - 



-22- 



1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 



r =0.6 

mode 



-10 12 3 
logD 



-15 



-16 



-17 



o -18 



-19 



-20 



-21 





1 ■ 1 ■ 1 I 1 ■ 1 ■ 1 1 1 1 1 






r =0.4 






mode 






• 






• 






• • 






• 






• • 






• « 












• • 












»• . 






• • 






... 












• 












•• 






• • 












* * • 






* * 






• 






• 










1 




_l 



1 2 3 4 5 
logD 



Fig. 5. — Appearance of L — D relation when increasing data dispersion (quantified with parameter 
r modc)- It can be inferred from the plots that larger scattering leads to the decrease and even loss 
of the correlation. 
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Fig. 6. — Performance of geometric mean regression for fitting power-law L — D relation L v = CD S . 
This type of regression can be a valid tool for obtaining L — D relation only for a limited range of 
data dispersion and it is closest to theoretical value 5 = —(3 + 2. 
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Fig. 7. — The S — D plots at 1 GHz for Balmer-dominated (empty circles) and oxygen rich SNRs 
(triangles) - thin line for the former and thick line for the later sample. The fit lines are obtained by 
orthogonal offsets. As expected, remnants evolving in a dense environment lie above those evolving 
in low density medium due to higher synchrotron emission. 
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Table 1. Shell SNRs with known distances* 



Surface brightness Distance Diameter 



Catalog name 


Other name 


(Wm^Hz-'sr- 1 ) 


(kpc) 


(pc) 


Referenc 


G4.5+6.8 a 


Kepler, SN1604, 3C358 


3.18e-19 


6.0 


5.2 


1 


Gil. 2-0.3 




1.93e-19 


4.4 


5.1 


7 


G18.8+0.3 a 


Kes 67 


2.66e-20 


12.0 


47.7 


2 


G21. 8-0.6 


Kes 69 


2.60e-20 


5.2 


30.3 


3 


G23.3-0.3 


W41 


1.45e-20 


4.2 


33.0 


2 


G27.4+0.4 


Kes 73, 4C-04.71 


5.64e-20 


8.65 


10.1 


2 


G31.9+0.0 a 


3C391 


1.03e-19 


8.5 


14.6 


2 


G33.6+0.1 a 


Kes 79 


3.31e-20 


7.0 


20.4 


4 


G41. 1-0.3 


3C397 


2.94e-19 


10.3 


10.0 


5 


G43.3-0.2 a 


W49B 


4.77e-19 


10.0 


10.1 


2 


G46.8-0.3 a 


HC30 


9.53e-21 


7.8 


33.7 


2 


G53.6-2.2 a 


3C400.2, NRAO 611 


1.30e-21 


2.8 


24.8 


2 


G54.4-0.3 a 


HC40 


2.63e-21 


3.3 


38.4 


6, 7 


G55. 0+0.3 




2.51e-22 


14.0 


70.5 


2 


G65. 1+0.6 




1.84e-22 


9.0 


175.6 


2 


G74.0-8.5 a 


Cygnus Loop 


8.59e-22 


0.54 


30.1 


8 


G78.2+2.1 a 


7 Cygni, DR4 


1.34e-20 


1.20 


20.9 


9 


G84.2-0.8 a 




5.17e-21 


4.50 


23.4 


7, 10 


G89.0+4.7 a 


HB21 


3.07e-21 


0.8 


24.2 


2 


G93.3+6.9 


DA 530, 4C(T)55.38.1 


2.51e-21 


2.2 


14.9 


7 


G93.7-0.2 


CTB 104A, DA 551 


1.53e-21 


1.5 


34.9 


2 


G94.0+1.0 


3C434.1 


2.61e-21 


5.2 


41.4 


2 


G96.0+2.0 




6.68e-23 


4.0 


30.3 


2 


G108.2-0.6 




3.19e-22 


3.2 


57.2 


2 


G109.1-1.0 


CTB 109 


4.22e-21 


3.2 


26.1 


11 


Gill. 7-2. l a 


Cassiopeia A, 3C461 


1.64e-17 


3.4 


4.9 


2 


G114.3+0.3 




1.67e-22 


0.7 


14.3 


2 


G116.5+l.l a 




3.14e-22 


1.6 


32.2 


2 


G116.9+0.2 a 


CTB 1 


1.04e-21 


1.6 


15.8 


2 


G119.5+10.2 a 


CTA 1 


6.69e-22 


1.4 


36.7 


2 


G120.1+1.4 a 


Tycho, 3C10, SN1572 


1.32e-19 


4.0 


9.3 


12 


G127.1+0.5 


R5 


8.92e-22 


1.25 


16.4 


2 
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Table 1 — Continued 







Surface brightness 


Distance 


Diameter 




Catalog name 


Other name 


(Wm^Hz - ^- 1 ) 

V / 


(kpc) 


(PC) 


Reference 


G132.7+1.3 a 


HB3 


1.06e-21 


2.2 


51.2 


2 


G156.2+5.7 a 




6.22e-23 


1.0 


32.0 


13 


G160.9+2.6 a 


HB9 


9.85e-22 


0.8 


30.2 


14 


G166.0+4.3 a 


VRO 42.05.01 


5.47e-22 


4.5 


57.4 


2 


G180.0-1.7 


S147 


3.02e-22 


0.62 


32.5 


2 


G189.1+3.0 


IC443, 3C157 


1.22E-20 


1.5 


19.6 


7 


G205.5+0.5 a 


Monoceros Nebula 


4.98e-22 


1.2 


76.8 


2 


G260.4-3.4 a 


Puppis A, MSH 08-44 


6.52e-21 


2.2 


35.1 


2 


G290.1-0.8 


MSH 11-61A 


2.38e-20 


7.0 


33.2 


2 


G292.2-0.5 




3.51e-21 


8.4 


42.3 


2 


G296.5+10.0 a 


PKS 1209-51/52 


1.23e-21 


2.1 


46.7 


15 


G296.8-0.3 


1156-62 


4.84e-21 


9.6 


46.7 


2 


G309.8+0.0 a 




5.39e-21 


3.6 


22.8 


7 


G315.4-2.3 a 


RCW 86, MSH 14-63 


4.18e-21 


2.3 


28.1 


2 


G327.4+0.4 


Kes 27 


1.02e-20 


4.85 


29.6 


2 


G327.6+14.6 a 


SN1006, PKS 1459-41 


3.18e-21 


2.2 


19.2 


2 


G332.4-0.4 a 


RCW 103 


4.21e-20 


3.1 


9.0 


2 


G337.0-0.1 


CTB 33 


1.34e-19 


11.0 


4.2 


2 


G337.8-0.1 


Kes 41 


5.02e-20 


11.0 


23.5 


2 


G338.3-0.0 




1.57e-20 


11.0 


25.6 


20 


G340.6+0.3 




2.09e-20 


15.0 


26.2 


2 


G344.7-0.1 




3.76e-21 


6.3 


18.3 


19 


G346.6-0.2 




1.88E-20 


11.0 


25.6 


16 


G348.5+0.1 a 


CTB 37A 


4.82e-20 


7.9 


34.5 


21 


G348.5+0.1 a 


CTB 37B 


1.35e-20 


13.2 


65.3 


21 


G349.7+0.2 a 




6.02e-19 


18.4 


12.0 


2 


G352.7-0.1 




1.25e-20 


7.5 


15.1 


4 


G359.1-0.5 a 




3.66e-21 


7.6 


53.1 


17, 18 



* Direct distance estimates, which can be inferred from proper motions, shock and radial ve- 
locities, HI absorption and polarization, association or interaction with HI, HII and CO molecular 
clouds, OB associations, pulsars, X-ray observations, optical extinction and low frequency radio 
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absorption. 

a SNRs also belonging to Case k Bhattacharya's (1998) calibration sample 

References. — (1) Chiotellis et al. 2012; (2) Green 2009; (3) Xin et al. 2009; (4) Giacani et al. 
2009; (5) Jiang et al. 2010; (6) Junkes et al. 1992; (7) Case k Bhattacharya 1998; (8) Blair k 
Sankrit 2005; (9) Uchiyama et al 2002; (10) Feldt k Green 1993; (11) Kothes k Foster 2012; (12) 
Asami Hayato et al. 2010; (13) Xu et al. 2007; (14) Leahy k Tian 2007; (15) Giacani et al. 2000; 
(16) Hewitt et al. 2009; (17) Uchida et al. 1992b; (18) Uchida et al. 1992a; (19) Giacani et al. 
2011; (20) Castelletti et al. 2011; (21) Tian k Leahy 2012; 
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Table 2. The £ — D fit parameters for calibration sample of 60 Galactic SNRs and different 
regression methods. We give f3 and A in relation £„ = AD~@ as well as their errors. 



Average 



Fit 


P 


A/3 


log A 


A log A 


fractional error 


OLS(F|A) 


2.2840 


0.3282 


-17.0028 


0.4895 


0.824 


OLS(A|Y) 


5.0375 


0.7465 


-13.1820 


1.0162 


0.469 


Orthogonal 


4.8161 


0.7218 


-13.4877 


0.9819 


0.472 


Bisector 


3.1741 


0.3335 


-15.7641 


0.4746 


0.575 


Arithmetic mean 


3.7079 


0.4282 


-15.0209 


0.5817 


0.522 


Geometric mean 


3.4178 


0.3536 


-15.4256 


0.4898 


0.547 
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Table 3. Distances to shell SNRs calculated from our £ — D relation 



Flux density Diameter Distance 





Ot.VlPT* TlflTTlP 


(Jv) 


(uc) 




GO.0+0.0 


Sgr A East 


100.0 


8.0 


9.3 


GO.3+0.0 




22.0 


18.4 


5.8 


Gl.0-0.1 




15.0 


17.5 


7.5 


Gl.4-0.1 




2.0 


28.8 


9.9 


Gl.9+0.3 




0.6 


17.1 




G3.7-0.2 




2.3 


30.5 


8.5 


G3.8+0.3 




3.0 


33.6 


6.4 


G4.2-3.5 




3.2 


39.7 


4.9 


G4.5+6.8 


Kepler, SN1604, 3C358 


19.0 


11.2 


12.9 (6.0)' 


G4.8+6.2 




3.0 


33.6 


6.4 


G5.2-2.6 




2.6 


34.6 


6.6 


G5.5+0.3 




5.5 


26.4 


6.8 


G5.9+3.1 




3.3 


34.4 


5.9 


G6.1+0.5 




4.5 


28.5 


6.7 


G6.4+4.0 




1.3 


49.6 


5.5 


G6.5-0.4 




27.0 


21.6 


4.1 


G7.0-0.1 




2.5 


32.4 


7.4 


G7.2+0.2 




2.8 


28.9 


8.3 


G7.7-3.7 


1814-24 


11.0 


28.1 


4.4 


G8.3-0.0 




1.2 


23.1 


17.7 


G8.7-5.0 




4.4 


36.1 


4.8 


G8.7-0.1 


W30 


80.0 


25.1 


1.9 


G8.9+0.4 




9.0 


30.3 


4.3 


G9.7-0.0 




3.7 


28.1 


7.5 


G9.8+0.6 




3.9 


27.1 


7.8 


G9.9-0.8 




6.7 


24.3 


7.0 


G10.5-0.0 




0.9 


27.5 


15.8 


Gll.0-0.0 




1.3 


31.3 


10.8 


Gil. 1-1.0 




5.8 


27.1 


6.3 


Gil. 1-0.7 




1.0 


31.4 


12.3 


Gil. 1+0.1 




2.3 


29.0 


9.1 



Table 3 — Continued 







Flux density 


Diameter 


Distance 


L/dldlOg IldlllC 


\J LIlcI Il<±IIl<J 


v J yj 






Gll.2-0.3 




22.0 


12.3 


10.5 (4.4) a 


Gll.4-0.1 




6.0 


21.1 


9.1 


Gil. 8-0.2 




0.7 


24.6 


21.1 


G12. 0-0.1 




3.5 


22.3 


10.9 


G12.2+0.3 




0.8 


27.2 


17.0 


G12.7-0.0 




0.8 


28.2 


16.1 


G13.5+0.2 




3.5 


18.6 


14.3 


G14.1-0.1 




0.5 


29.9 


18.7 


G14.3+0.1 




0.6 


26.5 


20.4 


G15.1-1.6 




5.5 


35.0 


4.5 


G15.4+0.1 




5.6 


27.2 


6.4 


G15.9+0.2 




5.0 


19.4 


11.3 


G16.0-0.5 




2.7 


29.4 


8.3 


G16.2-2.7 




2.0 


35.7 


7.2 


G16.4-0.5 




4.6 


27.0 


7.2 


G17.0-0.0 




0.5 


28.8 


19.8 


G17.4-2.3 




4.8 


34.3 


4.9 


G17.4-0.1 




0.4 


32.4 


18.6 


G17.8-2.6 




4.0 


35.6 


5.1 


G18.1-0.1 




4.6 


22.2 


9.6 


G18.6-0.2 




1.4 


25.2 


14.4 


G18.8+0.3 


Kes 67 


33.0 


18.6 


4.7 (12.0) a 


G19.1+0.2 




10.0 


31.1 


4.0 


G20.4+0.1 




3.1 


24.1 


10.3 


G21.0-0.4 




1.1 


29.6 


12.8 


G21.5-0.1 




0.4 


30.1 


20.7 


G21.8-0.6 


Kes 69 


69.0 


18.6 


3.2 (5.2) a 


G22. 7-0.2 




33.0 


24.0 


3.2 


G23.3-0.3 


W41 


70.0 


21.0 


2.7 (4.2) a 


G24.7-0.6 




8.0 


25.6 


5.9 


G27.4+0.0 


Kes 73, 4C-04.71 


6.0 


15.9 


13.7 (8.65) a 


G28.6-0.1 




3.0 


27.4 


8.7 
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Table 3 — Continued 







Flux density 


Diameter 


Distance 


Catalog name 


Other name 


(Jy) 


(pc) 


(kpc) 


G29 6-1-0 1 

VI Zi ,J . \J | \J . X 




1 5 

J. . '-> 


23.1 


15 9 


osn 7-i-i n 




fi 


31 


5 1 


osi 5— n r 




2 


3fi 5 


7 




3G3Q1 


24.0 


14.1 


8 2 f8 51° 


0^2 D— 4 Q 


3C39R 1 


22.0 


3fi fi 


2.1 


fiw 4-1-0 1 




2 

V/ . Z 


35 fi 


20.4 


0^2 8— fl 1 


Kps 78 


1 1 


25 3 

ZO . '.) 


5 1 


o^s 9— n fi 

V.' tJ ' ) . Z \ ) . \) 




3 5 


32 fi 


fi 2 




Kps 79 4G00 70 HG1 ^ 


99 

ZjZ . V 


1 7 7 

±1.1 


fi 1 (7 Q) a 


OSfi fi_(_2 fi 




7 


41 7 


9 7 






1 1 

X X . L/ 


98 1 


4 4 


Q41 i_n q 


3C397 


22.0 


11.4 


11 7 C10 3 s ) a 

XX. 1 1 IWiO 1 


r;42 S-un fi 




3 


37 8 


5 4 


G43 3-0 2 


W49B 


38.0 


10.4 


10 3 f 10 0) a 


Q43 qii k 




8 fi 

O . VJ 


44.2 


2 5 

Z. . 


riAK 7— n 4 




4.2 


34 1 


5 3 


VI i Vy . O VJ.O 


HC30 


14 


22.8 


5 3 (7 R) a 


Q40 o— n 7 


W51 


160.0 


18.5 


2.1 


G53 fi— 2 2 

VI • ) ' ) . \) Zj.Zj 


3C400 2 NR AO 61 1 

'.) V. x \ )\J • Z^ j ± > X IAV V 1 1 


8 

O . V/ 


34 1 


3 9 C2 8l a 


Q54 4—0 3 


HC40 

1 1 V t: W 


28 

ZjO . V/ 


29 fi 


2 5 f3 3) a 

Zj . tj 1 (J.U 1 


G55 O-t-D 3 




5 


47 5 


q 4 (44 nl* 1 


Q55 7+q 4 

V I ' ' J ■ 1 | ' • t: 




1.4 


43.3 


6.5 


G57 2-1-0 8 


4C21 53 


1.8 


31.6 


9.1 


Q50 5 in i 

VI ' 7 • / • 'J | V/ . 1 




3.0 


31.2 


7.2 


G65.1+0.O 




r r 

5.5 


50. 6 


z.6 (9.0J 


G65.3+5.7 




52.0 


56.6 


0.7 


G67.7+1.8 




1.0 


37.3 


9.5 


G69.7+1.0 




2.0 


33.9 


7.8 


G73.9+0.9 




9.0 


31.7 


4.0 


G74.0-8.5 


Cygnus Loop 


210.0 


37.1 


0.7 (0.54) a 


G78.2+2.1 


DR4, gamma Cygni SNR 


320.0 


21.3 


1.2 (1.2)° 


G82.2+5.3 


W63 


120.0 


29.0 


1.3 
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Table 3 — Continued 







Flux density 


Diameter 


Distance 


Catalog name 


Other name 


(Jy) 


(PC) 


(kpc) 


083 n_n S 

V.I (J ' ) . V I \J . 




1 

X . \J 


30 2 


1 3 1 

lO.l 


VJOt: . Zj V/.O 




1 1 


25 8 


5 (A 5"l a 




HB21 


220.0 


28.7 


9 (0 8) a 


V.I • / ' } . (J | VJ • > 7 


DA 530 40(1+55 38 1 


9.0 


30.3 


4.5 (2.2) a 




CTR 104A DA 551 


65 


33 


1 4 fl 5") a 


nqA n_i_i n 


SG434 1 


1 3 


2Q 6 


3 7 C5 2") a 


nqa n_i_9 n 




3 


62 


8 2 C4 0l a 


G108 2—0 6 




8.0 


45.3 


2 5 f3 2) a 


V.I 1 W . / • 1 X.U 


CTR 109 

V X XJ XL/iV 


22.0 


26 Q 

ZjL» . tJ 


3 3 f3 2l a 


Gill. 7— 2.1 


Cassioneia A 3C461 


2720.0 


5.1 


3.5 (3.4) a 


01 14 S-l-0 3 

V.I X 11. O ] . o 




5 5 


51 6 


2 5 (0 7l a 


01 1 6 5-1-1 1 

V.I 'J \ X . -L 




1 


45 4 


2 3 fl 6) a 


Ol 1 6 Q-l-D 2 

V.I is | VJ . Zj 


CTR 1 


8 


35 7 


3 6 fl 6") a 


r;i iq 5 -Lin 2 

V.I X X iV . <J ]^ 


CTA 1 


36 


39 


1 5 fl AY 


C12D 1-1-1 4 

V .I 1 Z.\ ) . 1 X . t: 


Tvrho 3C10 SN1572 


56 


1 3 4 

X (J.I 


5 8 f4 f]\ a 


CI 26 2+1 6 

V .I 1 Z.U. Zj^^ J- • 




6 


50 6 


2 5 

Zj .fj 


CI 27 1-t-O 5 

V.I X £j 1 . X | VJ • (J 


R5 


12.0 


36.8 


2 8 fl 25) a 


r;i 7+1 3 

V.I X (JZj . 1 . > ) 


HR3 


45 


35 5 


1 5 f2 2) a 

X . *J \ Zj . Zj 1 


CI 56 2-1-5 7 

V .I 1 • ) \ ) . Zj \ - J . 1 




5 


62 9 


2 fl 0l a 

Zj . V/ \ X . \J J 


0160 q_i_9 6 

V.I X \J\J . •Ji^jLj . \J 


HR9 


110 


36 1 

OU. X 


1 fO 8l a 


CI166 n_i_4 ^ 

V .I 1 V.) \ ) . \ ) t: . ' J 


VRO 42 05 01 


7 


40 6 


3 2 f4 5 N ) a 


r;i7q n_|_9 6 

V.I 1 1 U • \J | -LJ . \J 




7.0 


49.0 


2.4 


rnan n_i 7 

V.I X CJVJ . VJ X . ( 


S147 


65 


45 8 


9 fO 62 N l a 


HI 82 4-1-4 




1 2 

X . L. 


61 1 

«J X . X 


4 9 


G189. 1+3.0 


IC443, 3C157 


160.0 


21.8 


1.7 (1.5)° 


G192.8-1.1 


PKS 0607+17 


20.0 


41.4 


1.8 


G205.5+0.5 


Monoceros Nebula 


160.0 


41.4 


0.6 (1.2) a 


G206.9+2.3 


PKS 0646+06 


6.0 


43.8 


3.1 


G260.4-3.4 


Puppis A, MSH 08-44 


130.0 


24.6 


1.5 (2.2) a 


G261.9+5.5 




10.0 


34.3 


3.4 


G266.2-1.2 


RX J0852. 0-4622 


50.0 


41.0 


1.2 


G272.2-3.2 




0.4 


46.9 


10.7 
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Table 3 — Continued 







Flux density 


Diameter 


Distance 


Catalog name 


Other name 


(Jy) 


(PC) 


(kpc) 


G279 0+1 1 




30.0 


41.3 


1.5 


C284 3—1 8 


MSH 10—53 


11.0 


29.1 


4.2 


C28fi 5—12 




1.4 


33 8 


3 


0289 7—0 3 




6 2 

\J.Li 


27 6 

Li 1 • \J 


fi 


0290 1 — 8 


MSH 1 1 — fil A 


42.0 


1 9 


4 (7 Q) a 


H2Q2 2—0 5 




7 


97 9 

Li 1 .Cf 


5 5 (R A\ a 


("120(1 1 —0 5 




8 


34 1 


3 Q 


C206 5+10 


PKS 1200—51/52 


48.0 


34.4 


1.5 (2.1) a 


ri2Qfi 8—0 3 


1 1 ^fi_fi9 

X X rJKJ \jLi 


9 


96 2 


5 4 fQ fi"!" 


090,8 fi— 




5 


94 3 


8 


0299 2—2 9 




0.5 


43.7 


10.7 


0200 fi— 5 

v.l Li .7 ."7 . V/ L/.fJ 




1 

X . V7 


36 8 


7 


0301 4—1 o 




2 1 

Li . X 


43 9 


5 9 


0302 3-I-0 7 




5 


99 6 


fi 


0304 fi_i_n 1 


Kps 1 7 

IVCO X i 


14 


1 7 8 

X i .O 


7 fi 


G308 1 —0 7 




1.2 


35 5 


4 


G309 2—0 fi 

V.l tJ V >7 . Li L/.LJ 




7 


25.2 


fi 4 


0309 8-1-0 




1 7 


95 6 


4 f3 fi"l a 

t:.V7 IO.u ( 


0310 fi— o 3 


Kps 20R 


5 


91 9 

Li ± . iJ 


4 


Q31 n 8—0 4 


Kes 20A 


6.0 


24.8 


7.1 






3.0 


20.1 


13.8 


C312 4—0 4 




45.0 


26.3 


2.4 


031 2 5—3 




3 5 


33 3 


fi 


Q315 4_2 3 


RCW 86 MSH 14—63 


49.0 


26.9 


2 2 (2 3l a 


G315.9— 0.0 




0.8 


44.6 


8.2 


G316.3-0.0 


MSH 14-57 


20.0 


24.0 


4.1 


G317.3-0.2 




4.7 


25.2 


7.9 


G321.9-0.3 




13.0 


29.3 


3.8 


G323.5+0.1 




3.0 


29.5 


7.8 


G327.2-0.1 




0.4 


30.1 


20.7 


G327.4+0.4 


Kes 27 


30.0 


22.5 


3.7 (4.85) a 


G327.4+1.0 




1.9 


33.3 


8.2 
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Table 3 — Continued 



Catalog name 


Other name 


Flux density 

(Jy) 


Diameter 
(pc) 


Distance 
(kpc) 


G327 6+14 6 


SN1006 PKS 1459—41 


19.0 


28.5 


3 3 (2 2) a 


ns30 04-15 n 


T ill nnc T ifinn 

XJ IXJJ LXo UVJVJ^J 


350 

' ) -J \ ) . yj 


32 6 


6 

V4 . U 


0330 24-1 n 




5 

KJ . KJ 


24.9 


7 8 


(4332 D4-D 2 




8 


23.4 


6 7 


(4332 4— D 4 


RCW 103 

1 V V_ V V X V4 >_) 


28 

Zjkj . yj 


1 6 9 

±yj . u 


5 8 f 3 1 ) a 

KJ . O ( O . X ( 


(4332 44-fl 1 

V_4 >44J ^_ . t: 1 V4 . X 


MSH 16—51 Kes 32 

IVXkJXX X VJ (Jlj (J^j 


26 


20.2 


4 6 


(4332 5—5 6 




2 


47 7 


4 7 


G335 24-0 1 

V.I ' ) > } KJ . jLj 1 VJ • J. 




16.0 


25.5 


4.2 


(4336 74-0 5 




6 


24 7 


7 2 


(4337 0—0 1 


CTR 33 

v__ X 1 4 (JO 


1 5 


14.2 


32 5 C1 1 n"i a 


G337 2-0 7 




1.5 


24.8 


14.2 


(4337 34-1 


Kes 40 

11CD t:V4 


1 6 


21.3 


5 5 

fX . O 


cm 7 8—0 i 


Kps 41 


1 8 
xo . yj 


1 6 3 


7 6 (4 1 0l a 

i . \J i x x . yj j 


(4S38 i_i_n 4 




4 


29 5 

Zj<J .KJ 


6 8 

y> . o 


V4 ' ) >_) O . ' ) \J ,\J 




7 


20.4 


8 8 (4 1 0l a 

o . o i x x . yj j 


QWR 54-0 1 




12.0 


19.2 


7 3 


(4^4(1 4-1-0 4 




5 


22 3 


9 1 

t/. X 


V_.I »X— tVJ . Un^VJ . O 




5 

KJ . V4 


19 5 


11 2 f15 f)1 a 

xx. £j \±Kj.yj) 


(4^41 q_n q 

V_.I (J41 . • 7 \J . ' ) 




2 5 

. KJ 


23 8 


1 1 7 

X X . 1 


C342 0-0 2 




3.5 


26.1 


8.6 


C342 14-0 9 




0.5 


37.3 


13.5 


(4343 i _n 7 




7 8 


31 

KJ X . VJ 


4 5 

t: . KJ 


G344 7-0 1 




2.5 


27.5 


9 5 f6 3) a 

kJ 'KJ \ \J • KJ J 


(4345 7—0 2 




6 


29 9 


1 7 1 

X 1 . X 


G346.6— 0.2 




o n 

8.0 


19.9 


o r /i i n\i7, 

8.5 (11. 0) 


G348.5-0.0 




10.0 


20.8 


7.2 


G348.54-0.1 


CTB 37A 


72.0 


16.5 


3.8 (7.9) a 


G348.7+0.3 


CTB 37B 


26.0 


21.3 


4.3 (13.2) a 


G349. 2-0.1 




1.4 


27.3 


12.8 


G349.7+0.2 




20.0 


9.9 


15.2 (18.4) a 


G350.0-2.0 




26.0 


31.5 


2.4 


G351.7+0.8 




10.0 


25.1 


5.4 



Table 3 — Continued 



Flux density Diameter Distance 



Catalog name Other name 


(Jy) 


(PC) 


(kpc) 


G351.9-0.9 


1.8 


29.9 


9.9 


G352.7-0.1 


4.0 


21.6 


10.7 (7.5) a 


G353.6-0.7 


2.5 


42.9 


4.9 


G353.9-2.0 


1.0 


36.8 


9.7 


G354.8-0.8 


2.8 


34.8 


6.3 


G355.4+0.7 


5.0 


34.6 


4.8 


G355.6-0.0 


3.0 


22.9 


11.3 


G355.9-2.5 


8.0 


24.2 


6.4 


G356.2+4.5 


4.0 


36.2 


5.0 


G356.3-0.3 


3.0 


25.2 


9.9 


G356.3-1.5 


3.0 


33.1 


6.6 


G357.7+0.3 


10.0 


29.6 


4.2 


G358.0+3.8 


1.5 


52.3 


4.7 


G358.1+0.1 


2.0 


38.1 


6.5 


G358.5-0.9 


4.0 


31.0 


6.3 


G359.0-0.9 


23.0 


24.6 


3.7 


G359.1-0.5 


14.0 


27.7 


4.0 (7.6) a 



* Using our £ — D relation for this SNR makes no sense, since it is 
reliably established that Gl.9+0.3 is increasing in its flux density. 

a SNRs belonging to our new calibration sample from Table 1. Val- 
ues in brackets represent directly obtained distances. 
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Table 4. Basic properties of the 28 Galactic SNRs evolving in dense ISM i.e. associated with 

molecular clouds 









Surfcicc brightness 


DicLinct ci' 


l. ci 1" p 1 t~i <y riamr 1 


( n Yl +- 1 c iTl CI 1 VI CI TY1 O 


1 i V 1LIL11LL 


fWm" 2 H7~ 1 sr" 1 'l 

l VV 111 Xlzj d! } 


(nc) 


G18.8+0.3 a 


Kes 67 


CO MA & LB, CO ratio 


2.66c-20 


47.7 


G21. 8-0.6 


Kcs 69 


OH, CO MA & LB, HCO+, H 2 


2.60e-20 


30.3 


G23.3-0.3 


W41 


HI MA & CO RC, extended TeV 


1.45e-20 


33.0 


G31.9+0.0 a 


3C391 


OH, molecular MA & LB (CO, HCO+,CS), H 2 , NIR 


1.03e-19 


14.6 


G33.6+0.1 a 


Kcs 79 


CO MA, HCO+ MA, broad OH absorption 


3.31c-20 


20.4 


G41. 1-0.3 


3C397 


CO MA & LB 


2.94c-19 


10.0 


G43.3-0.2 a 


W49B 


H 2 MA, CO ratio 


4.77c-19 


10.1 


G54.4-0.3 a 


HC40 


CO MA & LB, IR MA 


2.63c-21 


38.4 


G78.2+2.1 a 


7 Cygni, DR4 


CO MA 


1.34c-20 


20.9 


G84.2-0.8 a 




CO MA 


5.17c-21 


23.4 


G89.0+4.7 a 


HB21 


CO MA & LB, CO ratio, H 2 , NIR 


3.07c-21 


24.2 


G94.0+1.0 


3C434.1 


CO RC 


2.61e-21 


41.4 


G109. 1-1.0 


CTB 109 


CO MA & LB 


4.22c-21 


26.1 


G111.7-2.1 a 


Cassiopeia A, 3C461 


H 2 CO absorption, IR RC, CO RC 


1.64c-17 


4.9 


G132.7+1.3 a 


HB3 


CO MA 


1.06c-21 


51.2 


G166.0+4.3 a 


VRO 42.05.01 


unusual shape, CO RC 


5.47c-22 


57.4 


G189. 1+3.0 


IC443, 3C157 


OH, CO ratio, H 2 , molecular MA & LB 


1.22E-20 


19.6 


G205.5+0.5 a 


Monoceros Nebula 


CO RC 


4.98c-22 


76.8 


G260.4-3.4 a 


Puppis A, MSH 08-44 


CO RC, OH(ncgativc) 


6.52c-21 


35.1 


G290.1-0.8 


MSH 11-61A 


CO RC 


2.38c-20 


33.2 


G332.4-0.4 a 


RCW 103 


IR MA & colors, NIR, H 2 & HCO+ MA 


4.21c-20 


9.0 


G337.0-0.1 


CTB 33 


OH 


1.34c-19 


4.2 


G337.8-0.1 


Kcs 41 


OH 


5.02c-20 


23.5 


G344.7-0.1 




IR RC & colors 


3.76c-21 


18.3 


G346. 6-0.2 




OH, H 2 , IR colors 


1.88E-20 


25.6 


G348.5+0.1 a 


CTB 37A 


OH, H 2 , IR MA 


4.82c-20 


34.5 


G349.7+0.2 a 




OH, CO MA & LB, CO ratio, H 2 , IR MA 


6.02c-19 


12.0 


G359.1-0.5 a 




OH, CO & H 2 MA, HCO+ & CS absorption 


3.66c-21 


53.1 



Evidence: chief evidence that suggests the interaction between SNR and molecular clouds. LB=line broadening, MA=morphology 
agreement, H 2 =vibrational/rotational lines of molecular hydrogen [e.g. H 2 1—0 S(l) line (2.12 um), H 2 0—0 S(0)— S(7) lines], 
NIR=Near-Infrared (e.g. [Fe II] line), OH=1720 MHz OH mascr, RC=rough morphological correspondence, etc. (taken from the List 
of Galactic SNRs Interacting with Molecular Clouds by Bing Jiang, bjiang@nju.edu.cn, 2010) 




kpc 



